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Abstract 

We consider the local chirality of near-zero eigenvectors from Wilson- 
Dirac and clover improved Wilson-Dirac lattice operators as proposed 
recently by Horvath et al. We studied finer lattices and repaired for 
the loss of orthogonality due to the non-normality of the Wilson-Dirac 
matrix. As a result we do see a clear double peak structure on lattices 
with resolutions higher than 0.1 fm. We found that the lattice artifacts 
can be considerably reduced by exploiting the biorthogonal system of 
left and right eigenvectors. We conclude that the dominance of in- 
stantons on topological charge fluctuations is not ruled out by local 
chirality measurements. 
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1 Introduction 

It is generally accepted that instanton field configurations play an important 
part in hadron phenomenology jjj]. The space-time structure of instantons 
has been the target of a variety of lattice investigations in the past 0. Re- 
cently Horvath et al. []|] have added a very interesting contribution to this 
discussion about the role of instantons. They studied the local chiral orien- 
tation of near-zero eigenfunctions on an ensemble of (rather coarse grained) 
quenched QCD vacuum configurations, but saw no evidence for instanton 
driven topological charge fluctuations. On the other hand, by applying the 
techniques of Horvath et al. to overlap Q instead of standard Wilson fer- 
mions, DeGrand and Hasenfratz arrived at the opposite conclusion ||. 
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Wilson-type actions have the merit of being ultralocal which is desirable 
for the study of local properties but at the cost of chiral symmetry breaking. 
On the other hand the overlap formulation has a remnant of chiral symmetry 
H but might wipe out local structures, as it implies couplings beyond nearest 
neighbours n, M. 

In view of the implications of the above findings for our understanding 
of the QCD vacuum structure we shall focus in this letter on possible lattice 
artifacts due to the use of Wilson-type operators. At finite lattice spacing a 
these operators are non-normal. In Section ||| we shall remind the reader of 
some peculiarities encountered in the eigenmode expansion due to this non- 
normality H . The latter induces a biorthogonal system of distinct right 
and left eigenvectors. This suggests to modify the local chiral orientation 
parameter X proposed in Ref. H. The resulting lattice improved form, 
X imp , is based on projecting the chiral reduction of the right eigenvectors 
onto the related left eigenvectors. In Section ||| we shall demonstrate this 
improvement for the case of QED2. In Section || we shall apply A imp to the 
case of QCD. 



2 Local Chirality in Non-normal Scenario 

Two commonly used discretizations of the Dirac operator are the Wilson- 
Dirac matrix, Dw, and the clover improved Wilson-Dirac matrix, Dew- 
Both matrices are non-normal 

rt/o. (1) 

This fact implies that both matrices are not diagonalizable by unitary trans- 
formations. Yet their diagonalization can be achieved by similarity trans- 
formations with non-unitary matrices 

A = S~ l DS , (2) 

with A being the diagonal matrix. The columns of S are then the right 
eigenvectors, \Ri), and the rows of S^ 1 the left eigenvectors, of D (see 
e.g. 0), 

D\Ri) = \i\Ri), 

(Li\D = A; (Li | . (3) 

It follows directly from S~ 1 S = I that the left and right eigenvectors are 
biorthogonal 

(L l \R J )=8 lj . (4) 

In contrast, the right (or left) eigenvectors by themselves do not form an 
orthogonal system. 
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This biorthogonality has to be taken into account when expressing an 
abstract operator A in terms of a matrix representation by use of eigenfunc- 
tions 

A = J2\R i ){L i \A\R j ){L j \. (5) 
»,j 

In the following we will consider the operator (1+75), which can be expanded 
accordingly 

l± l5 = J2\R i )(L i \l± l5 \R j )(L j \ . (6) 
hi 

(Actually, we will need only the diagonal elements =L 7s|-Rj)). Only in 
the continuum limit, where D becomes normal, the operators in Eqs. (||) 
and © can be expanded in terms of right eigenvectors only 

Fortunately, though, one encounters no duplication of the costs to com- 
pute the left eigenvectors of D\y and Dew because these operators obey the 
75-hermiticity property 

L>t = 75^75 • (7) 

This symmetry allows for a straightforward construction of the left eigen- 
vectors from the right ones: 

(R i \j 5 D = X*(R i \ l5 . (8) 

Hence the left eigenvector to the eigenvalue A* is nothing but the hermitian 
conjugate of the right eigenvector to Aj multiplied by 75 from the right: 

(L l \D = (R i \ l5 D = X*(L i \. (9) 

Horvath et al. method. In order to quantify the role of instantons in 
the topological charge fluctuations in QCD Horvath et al. define the local 
"chiral orientation" parameter X of an eigenmode \rp) at every lattice point 
(cf. Eq. (19) in 0), 

The absolute values \ip-\ and stand for [] 

h/>±| 2 = <^±h/'±) = (V'|i±75h/>>, (ii) 

where the scalar products are meant in colour and spin space, on an in- 
dividual lattice site. The analysis is carried out on 1 % of the lattice 
sites with maximal values (ip\ifi). As a result they produce histograms of 
X-distributions from a number of near zero-eigenmodes and gauge con- 
figurations. This is motivated by a continuum consideration according to 

1 We use ip- and ip+ instead of their ipL and ipR to avoid confusion with the left and 
right eigenvectors. 
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which QCD instantons produce distributions in the quantity X that peak 
at X = ±1. 

This procedure can be applied directly to normal Ginsparg- Wilson fer- 



mions where the choice of the bra- vector in Eq. (11) is not to be ques- 
tioned since the left and right eigenvectors coincide. The situation is less 
clear once the normality is lost. In our notation Horvath et al. start from 
right eigenmodes and usef] 

tan( i (l + X)j = |W+7s|W/2 . (12) 

For non-normal operators, however, the question arises whether it would 
be more economical to adjust to the above described biorthogonal structure 
of operator representations. This suggests to replace in Eq. (ITTj) 



(^±hkfc)-»(£i|l±T5|fli>, (13) 

for Wilson-type fermions. 

We emphasize that this modified prescription will coincide with the pre- 
vious one in the continuum limit, a — > 0, where the Dirac operator is normal. 
Nevertheless one might hope that a lattice adapted definition of X will sup- 
press lattice effects. With this motivation we propose here to make use of 
an "improved" quantity X imp , 

tan^(l + X imp )j = |(Lj|1 + 75| ^ )|1/2 . (14) 

The selection of test points is consistently chosen with respect to \(Li\Ri)\ 
(instead of (Ri\Ri)). 

Throughout this letter, we shall refer to the two definitions as follows: 

{Ri\l ± j 5 \Ri) ^ ^definition, (15) 
(Li\l ± j 5 \Ri) 9* Li? definition. (16) 



3 Results for QED2 

Let us consider two-dimensional quantum electrodynamics (QED2) in order 
to probe the concept of the "improved" local chiral orientation parameter 
X imp . In QED2 the gauge field is (anti)self-dual and one should therefore 
observe peaks in the chiral orientation near plus and minus one. 

In order to test the LR definition, Eq. (|l~4j), we perform an analysis on 
a 24 2 lattice for the /3-values 6.0, 2.0 and 1.0. We determined the complex 
eigenmodes in the physical branch with |ImA| < 0.3. The resulting his- 
tograms for 500 quenched gauge field configurations for each (3 are shown in 

2 I. Horvath, private communication 



4 



/3 = 6, RR-def inition 



3=6, LR-def inition 



5000 
4000 
3000 
2000 



5000 
4000 
3000 
2000 
1000 



6000 
5000 
4000 
3000 
2000 
1000 



-1 -0.5 0.5 
13 = 2, RR-def inition 



= 1, RR-def inition 



5000 
4000 
3000 
2000 
1000 



5000 
4000 
3000 
2000 
1000 



6000 
5000 
4000 
3000 
2000 
1000 



-1 -0.5 0.5 1 

13 = 2, LR-def inition 



13 = 1, LR-def inition 



-0.5 



0.5 



-1 -0.5 0.5 1 



Figure 1: Histograms of the local chiral orientation parameter for (complex) 
near-zero eigenvectors in QED2 for standard Wilson-Dirac fermions. The 
two definitions, RR and LR, have been used at (3 = 6.0, 2.0 and 1.0. 

the right column of Fig. |l|. Throughout our /3-range we are able to observe 
a very clear double peak structure. For comparison, we also included the 



histograms based on the RR definition, Eq. (12), for the same sets of con- 
figurations. We find that the signal becomes almost invisible at j3 = 2.0 and 
vanishes below. We emphasize that the results from both definitions tend 
to coincide for small lattice spacings, as they should. On the other hand, 
one can easily be misled by using the RR definition at too coarse lattice 
spacings. 

The reason behind this remarkable difference between the two schemes 
can be traced to the loss of orthonormality of the right eigenvectors. For 
if the left and hermitian conjugate right eigenmodes coincide, the following 
"cosine" 

{Li Ri} = v^jj^j) (i7) 

V (Li\Li)y/ (Ri\Ri) 
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Figure 2: Needle plots of the quantity {Li,Ri} (Eq. (Jl?])) computed for 
entire spectra on 16 2 lattices at j3 = 6.0 (left) and 1.0 (right). The lengths 
of the vertical needles represent the values of {Lj, R{\ for eigenmode i, while 
their positions are given in the complex Aj plane. 
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Figure 3: Same as Fig. |l] but with clover improved Wilson-Dirac fermions. 



would be equal to one. In Fig. ^ we show this quantity in form of needle plots 
for the entire spectrum of eigenmodes for a 16 2 lattice at j3 = 6.0 (left) and 
1.0 (right) in form of vertical lines positioned at the complex eigenvalues. 
For a normal operator, {Li, Ri} is identical to 1 (for all i). Hence deviations 
from one give us a measure for non-normality. 

The values of {Li, Ri} are found to be significantly closer to one for the 
larger /3-value. This holds in particular in the physical branch (near the left 
edge of the box). In fact, at the lower (5- value, we are too far away from one 
as to still rely on the orthogonality of the right eigenvectors alone. 

As another check of the consistency of this picture we have examined 
to what extent an improved action helps to enhance the signal. To this 
end we have computed the chiral orientation for the same set of quenched 
configurations but with clover improved Wilson-Dirac fermions. The results 
are gathered in Fig. || which is to be compared directly to Fig. p]. We 
find that the improvement on the signal is fully consistent with the idea of 
accelerated convergence (Ri\ — > (Li\. 



4 Results for QCD 

Encouraged by our results in QED2 we shall now turn to quenched QCD. 
To compare with the results of Ref. || , we computed 50 near-zero (complex) 
eigenmodes on 30 configurations of size 12 3 x 24 at the same /3-value, (3 = 
5.7. The resulting chirality distributions are plotted in Fig. ||, for both 
the RR and the LR definitions. While the -R-R-result shows only a plateau 
very similar to the one seen in Ref. Q, we do find a slight indication of a 
double peak structure in the LR case. As a next step, we proceed to finer 
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Figure 4: The histograms for X for the two definitions at f3 = 5.7 for QCD 
for near-zero eigenvectors of the Wilson-Dirac matrix, Dw ■ 

lattices, corresponding to f3 = 5.8, 6.0 and 6.1, on 16 3 x 32 configurations. 
The results from the RR and LR definitions are confronted in Fig. ||. We 
see unambiguous evidence for a deepening of the valley between two clear 
peaks, which appear to move outwards, towards ±1. This is confirmed by 
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the results from the RR definition, though with quite some delay when 
increasing j3. In order to achieve the same signal quality one has to work 
with = 6.1 (RR), instead of (3 = 5.8 {LR). 

To complete the study we repeated the computation with clover improve- 
ment, at /3 = 6.0. The histograms are shown in Fig. ^. As in the QED2 case 
we observe a markedly better separation of the peaks in the RR situation. 
There is also some further signal improvement in the LR case. 
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Figure 5: The histograms for X for the two definitions at f3 = 5.8, (3 = 6.0 
and (5 = 6.1 for QCD for near-zero eigenvectors of the Wilson-Dirac matrix, 
Dw 



5 Summary and conclusion 

In their debate about the conclusiveness of the local chirality data presented 
in Ref. ||], the authors of || have suspected that the artifacts and the use 
of a non-chiral lattice action might hide the signal for correlations between 
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Figure 6: The histograms for X at (3 = 6.0 in QCD for the two defini- 
tions using the near-zero eigenvectors of the clover-improved Wilson-Dirac 
operator, Dew- 



instanton positions and wave function localizations. This is precisely what 
we confirm in our analysis. 

To clarify the situation, we studied finer lattices and repaired for the loss 
of orthogonality due to the non-normality of the Wilson-Dirac matrix. As 
a result we do see a clear double peak structure on lattices with resolutions 
higher than 0.1 fm {(5 > 6.0). We found that the lattice artifacts can be 
considerably reduced by exploiting the biorthogonal system of left and right 
eigenvectors. 

We conclude that the dominance of instantons on topological charge 
fluctuations is not at all ruled out by local chirality measurements. 

In order to consolidate this picture it appears worthwhile to investigate 
which fermion formulation is actually the most appropriate in this context: 
although the overlap fermions seem to yield convincing peaks it is not clear 
whether less local actions might wipe out the fluctuations that one is after. 
But now the door is open for investigations of the local chiral structure, 
relying on ultralocal probes. This will help to elucidate the significance of 
the overlap signal. 
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